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ABSTRACT 

To explain the observed anomalies in stellar populations within globular clusters, many 
globular cluster formation theories require two independent episodes of star formation. A fun¬ 
damental prediction of these models is that the clusters must accumulate large gas reservoirs 
as the raw material to form the second stellar generation. We show that young clusters con¬ 
taining the required gas reservoir should exhibit the following observational signatures: (i) 
a dip in the measured luminosity profile or an increase in measured reddening towards the 
cluster centre, with Ay > 10 mag within a radius of a few pc; (ii) bright (sub)mm emission 
from dust grains; (iii) bright molecular line emission once the gas is dense enough to begin 
forming stars. Unless the IMF is anomalously skewed towards low-mass stars, the clusters 
should also show obvious signs of star formation via optical emission lines (e.g. Ha) after 
the stars have formed. These observational signatures should be readily observable towards 
any compact clusters (radii of a few pc) in the nearby Universe with masses > 10 6 M ( . ; and 
ages <100 Myr. This provides a straightforward way to directly test globular cluster formation 
models which predict large gas reservoirs are required to form the second stellar generation. 
The fact that no such observational evidence exists calls into question whether such a mecha¬ 
nism happens regularly for YMCs in galaxies within a few tens of Mpc. 
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1 INTRODUCTION 


Globular clusters (GCs) were once thought to contain stellar pop¬ 
ulations with the same age, metallicity and chemical abundances. 
However, it has become clear over the last few decades that this pic¬ 
ture is over simplified. When observed with sufficient photometric 
and spectroscopic precision, all GCs display anomalies, from mul¬ 
tiple discrete sequences in colour-magnitude di agram s to variations 
in the abundances of light elements te.g. lGratton et al.l2012li . These 
anomalies offer an opportunity to constrain the formation history of 
GCs - a challenging task given their formation at early epochs of 
the Universe. 

Several GC formation models have been put forw ard to ex¬ 
plain the origin of the observed anomalies (e.g. rDecressinetal.1 
20071: ID’Ercole et alj 120081: Ide Mink et all 120091: IConrovI 1201 ± 
Bastian et alJl2013h . Most of these invoke two or more episodes, 
or ‘generations’, of star formation, in which the second generation 
of stars forms from gas that has been enriched by material from 
certain subsets of stars in the first generation. A fundamental pre- 
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diction of the 'AGB' model scenarios (in which the polluters are 
asymptotic giant branch (AGB) stars) is that GCs must accumulate 
a large gas reservoir after the first generation have finished forming, 
in order to have sufficient material to form the second gene ration 
of stars dD’Ercole et al .1120081 l20ld : IConrov & Snereefcoi ill . 

In this paper we attempt to constrain the expected properties of 
this gas reservoir. We then consider how the gas may affect the ob¬ 
served properties of young GCs in the phase between the formation 
of the first and later generation of stars. As GC formation models 
make no distinction between the physics of star/cluster formation 
at the present day and earlier epochs of the Universe (when GCs 
were forming), these same effects should be seen towards clusters 
of similar mass, density and age in the local Universe. 

2 EXPECTED PROPERTIES OF GAS DESTINED TO 
FORM THE SECOND GENERATION OF STARS 

2.1 An idealised GC 

We first consider an idealised GC which forms two distinct stellar 
populations. The first generation of stars has a total initial mass of 
Ml, that all form at a time t\ = 0, i.e. as a single stellar population 
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(SSP). These stars are contained within a radius, R*. At some time 
later, f| tart , the second generation of stars begin to form. Star for¬ 
mation proceeds for a time interval of Af 2 , until the second gener¬ 
ation of stars finishes forming at t| nd = f| tart + Af 2 . The resulting 
mass of the second generation of stars, M|, is a fraction, F 2 , of the 
first generation (i.e. M| = T 2 M*). The time-averaged star forma¬ 
tion rate required to form the second generation of stars, SFRf v is 
then given by SFR| V = F 2 M* /A£ 2 - R 2 is defined as the radius 
containing the second generation of stars at 1 2 nd . In order to form a 
mass M| of stars, a total gas reservoir of mass, Mj otGas = M| It 
is required, where t is the star formation efficiency. 

To match the observed anomalies, the two-generation GC for¬ 
mation models have two distinct gas components - a component of 
enriched material expelled from stars within the cluster, and a com¬ 
ponent of ‘pristine’ (i.e. unprocessed) material. In AGB models, 
the pristine gas must be accreted from outside the cluster bound¬ 
ary. The mass and spatial distribution of gas within the GC after 
the first generation of stars have formed must therefore be time 
variable fsee lD’Ercole et af]|20ld ). In order to match the observed 
light element abundance variations, the two gas components must 
be at least partially mixed. This implies that both gas components 
have a high volume-filling factor, i.e. the gas distribution is fairly 
uniform. We discuss the implications of relaxing this assumption in 
Sj3] In the absence of numerical modelling to study the gas accumu¬ 
lation and star formation process in detaiQ, it seems reasonable to 
assume that the volume containing the gas before star formation be¬ 
gins (i.e. at time f| tart ) is similar to the final volume containing the 
second generation starfl Following these considerations, we make 
the simplifying assumptions that: (i) the radius containing the gas, 
R gas , is equal to R£; (ii) R gas remains fixed from f! tart to f| nd ; and 
(iii) the gas within this radius is uniformly distributed. 

The instantaneous gas mass within the cluster, M 2 nstGas , an[ j 
the insta ntaneous star formation rate , SFRj; nst , will vary with t ime 
(see e.g. iD’Ercole et alJl2oToL 1201 li t IConrov & Spergelll201 lh . In 
order to calculate the expected range in jVfi nstGas an d SFR 2 nst over 
the formation period of the second generation of stars, we consider 
two extreme formation scenarios. In the ‘continuous’ scenario, the 
gas accumulation rate is constant between f| tart a nd t | nd . The in¬ 
stantaneous mass accretion rate, M, is then M = Mj otGaa /At = 
M 2 /eAt. The instantaneous star formation rate is constant at a 
value of SFR 2 nst = tM. In the second ‘burst’ scenario, the full 
gas reservoir, Mj otGas , is in place at f| tart . This gas is then all 
converted into stars with a star formation efficiency, t, in a single 
star formation event lasting a time, Af 2 - We adopt a high value 
of t = 0.5 to provide a conservative lower limit to the required 
gas reservoir (and hence a lower limit to the extinction and mm- 
continuum flux density later in the paper). 

Having constructed this template, we now consider appropri¬ 
ate values for the parameters defined above. Due to stellar and 
dynamical evolution over close to a Hubble time, the initial mass 
(M*) and radius (Rj) of the first stellar generation in GCs are not 
empirically well constrained. Indeed, inferring M| and F 2 from the 
present-day stellar populations in GCs is mo del d ependent. For ex¬ 
ample, to solve the ‘mass budget’ problem, iD’Ercole et ali ( t2008i) 
require GCs to have been a factor ~ 10 more massive at birth than 
observed today. However, the mass of stars formed in the second 


generation. M|, is constrained much more robustly. Present-day 
second generation stars are confined t o a smaller volume than that 
of present-day first generation stars iLardo et~aD |20 1 ll ). In order 
to minimise the mass budget problem, it is generally assumed that 
the majority of second generation stars have remained bound since 
forming. In other words, modulo stellar evolution, the measured 
present-day mass of second generation stars is equal to M| . It fol¬ 
lows that the present-day mass of second generation stars is also a 
robust way to estimate the required gas reservoir. 

Although the initial size of GCs is not well constrained, a sur¬ 
prisingly uniform, present-day, half-mass radii of 3 pc is observed 
for GCs both in the Milky Way (MW ) GCs and in external galax¬ 
ies dHarri 3 l 19961 : iMasters et all koiQ). As mentioned above, sec¬ 
ond generation stars are confined to smaller radii than first gener¬ 
ation stars . We therefore adopt an upper limit for R 2 of 2 pc (see 
lLardo et al J20 1 ll . for observed spatial distributions of first and sec¬ 
ond generation stars). 

Limits for gas accumulation timescales in AGB GC formation 
models are set by the window in which SNe II in the first generation 
have finished and the time at which SNe la from the first genera¬ 
tion begin. The justification for these limits are that SNe II clear 
out any remaining pristine gas that did not end up in stars when 
the first generation formed, and the SNe la can clear any gas that 
builds up in the GC before star formation can begin. This gives an 
upper limit on the star formatio n timescale of t| tart ^ 30 Myr to 
f| nd ~ 100 Myr ( ID’Ercole et al.l[2008lFl which we adopt for the 
‘continuous’ star formation scenario. For the ‘burst’ scenario, the 
upper time limit, f| nd ~ 100 Myr, is the same. However, in this 
scenario, while the gas can continue to accumulate to high density, 
some mechanism must st op it from forming stars. In the model of 
IConrov & Spergel 1201 ll) . the high Lyman-Werner photon density 
at cluster ages < 10 s yr photodissociates molecular hydrogen (H 2 ), 
thereby suppressing star formation. At any given time it is assumed 
that the cluster contains a gas reservoir with a mass ~10% of the 
stellar mass. This helps the cluster to sweep up pristine gas from 
the surrounding environment. After t ~ 10 s yr the Lyman-Werner 
photon density drops precipitously, allowing the gas to catastroph¬ 
ically cool and form stars. In the ‘burst’ scenario, we adopt At 2 
to be ~ 1 Myr, corresponding to an upper limit of the observed age 
spread in young massive clusters (cf. lLongmore et ahll2014l) . 


2.2 Infall model 

The idealised scenario above is useful to constrain the expected 
range of g as properties in young GCs. We use the results of 
D^Ercole_etal. (20081 to consider more realistic gas profiles. 
D’Ercole et al. <20081) constructed a model of gas infalling towards 
the centre of a GC. This model is ideal for our purposes, as the in¬ 
falling gas reservoir is the raw material for the second generation 
of stars. We extracted an approximate density, p, and temperature, 
T, profile for the fiduciaJ_model _in t heir p aper, when the GC is 
at an age of 100Myr dD’Ercole et all 2008, Figure 1). From this 
figure we took the density and temperature structure to be power 
laws of the following forrrQ p(r) = 851 (r/pc) -2 ' 1 2 cm -3 and 


1 We consider more realistic gas profiles in § 12.21 

2 Observations show that young massive cluster progenitor clouds in the 
Milky Way are only factors of a few larger than the final stellar populations 
iLongmore et all20l4 Walker et al. 2015) 


3 We note that IConrov & Spergej |201 ll) adopt a longer (but not dir ectly 
specified) limit for f| nd of a few 100 Myr 

4 We tried several different parametric fits and the results are not sensitive 
to the exact form we chose. We therefore opted for the simplest, power law 
representation. 
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T(r) = 1548 (r/pc) 0 ' 34 K. These reproduce the density and tem¬ 
perature profiles over the range 0.04 pc < r < 60 pc. The lower 
limit to the fit_of 0.04 pc is set by the minimum radius in Figure 1 
of lD’Ercole et al.l J2008I ). We then constructed a spherical model of 
the gas from these profiles and used this to calculate the expected 
observational effect of this gas. To avoid unphysical temperatures 
and densities as r —> 0, we assumed the temperature and density 
at radii less than 0.04 pc to equal those at r = 0.04 pc - i.e. a flat 
temperature and density profile for r < 0.04 pc. In practice, this 
provides a lower limit to the expected extinction, gas emission and 
dust emission at r < 0.04 pc. 


3 IMPLICATIONS FOR THE OBSERVED PROPERTIES 
OF YOUNG GCS 



We now consider how the gas mass accumulation and subsequent 
star formation may affect the observable properties of young GCs. 
As our empirical understanding of GCs is gleaned almost ex¬ 
clusively from observations at optical and infrared wavelengths, 
we start by considering how gas accumulation and star forma¬ 
tion would affect such observations. One potential effect of the gas 
would be extinction and reddening of the light due to the dust as¬ 
sociated with the gas. With an estimate of the gas surface density 
from § [2] we can simply use a gas-to-dust ratio and an extinction 
law to calculate the expected extinction and reddening. To con¬ 
vert from the atomic gas column density, Nh, to V-band (wave¬ 
length of 5470 A) dust extinction, Ay, we start with a conversion of 
Av/Nh = 5.3 x 1 0~ 22 mag cm 2 H _1 , appropriate for Milky Way 
dust dDrain l20lH) . Making the common assumption that the gas- 
to-dust ratio scales linearly with metallicity (a reasonable assump¬ 
tion for different star formation environments over several orders of 
magnitude in metallicity: see lFisher et al.l20l4h . we can extrapolate 
the expected extinction to lower metallicity environments. 

FigureQ]shows the resulting expected extinction in the V-band 
as a function of the mass in the second generation of stars (i.e. 
as defined above) in the idealised GC (§ 12.11 - i.e. when the gas is 
confined to a radius of 2 pc). The black, red and blue lines represent 
gas-to-dust ratios of 10 2 , 10 3 and 10 4 , respectively. Assuming the 
gas-to-dust ratio scales linearly with metallicity, the lines represent 
metallicities of 1, 0.1, and 0.01 times solar. The solid and dashed 
lines shows the expected extinction for the ‘continuous’ and ‘burst’ 
scenarios, respectively. For a given gas-to-dust ratio (or metallicity) 
and MT the values between the solid and dashed lines show the 
range of expected extinction. It is clear that for gas-to-dust ratios 
^ 10 3 (i.e. ZIZq ^ 0.1) the V-band extinction is always >1 mag. 
As such, the gas and dust would be expected to have a significant 
effect on the observed optical magnitudes and colours of young 
GCs. 

To determine what the extinction profile might look like, we 
first projected the 3D volume density profile of the model con¬ 
structed in § 12.21 on to a 2D surface. We then converted the gas 
column density to a V-band extinction as a function of position us¬ 
ing the same conversion as above. The left panel of Figure[2]shqws 
the extinction map for the inner 30 pc of the lD’Ercole et alj 120081) 
model. The outer black contour shows that at radii <25 pc the pre¬ 
dicted extinction is >0.1 mag, so should be readily detectable in 
optical/infrared colour-colour and colour-magnitude diagrams. The 
right panel of Figure[2]zooms in to radii of <8 pc. At radii ~5 pc the 
extinction reaches 0.5 mag, and rises to >10 mag at pc-scales.This 
extinction would have a dramatic affect on optical/infrared obser¬ 
vations, causing a pronounced dip in the luminosity profile at small 


Figure 1. Expected V-band extinction as a function of mass in the second 
generation of stars. A/). The black, red and blue lines represent gas-to-dust 
ratios of 10 2 , 10 3 and 10 4 , respectively. Assuming the gas-to-dust ratio 
scales linearly with metallicity, the lines represent metallicities of 1,0.1, and 
0.01 times solar. The solid and dashed lines shows the expected extinction 
for the ‘continuous' and ‘burst’ scenarios, respectively 


radii. Such a dip should be easily detectable in spatially-resolved 
observations of young clusters. With multiple filter, high precision, 
spatially-resolved, optical/infrared photometry, the extinction may 
be detectable as increased reddening and extreme extinction to¬ 
wards the cluster centre. In summary, it should be possible to iden¬ 
tify any resolved, massive (> 10 6 Mq), compact (radii of a few pc 
and larger) clusters with large central gas reservoirsQ 

The dust causing the extinction will also be observable at 
longer wavelengths in emission. Figure [3] shows the expected flux 
density at a representative wavelength of 1 mm as a function of 
distance for gas of dif ferent masses, tempera tures and metallicities 
(calculated following Kauffmarm_etal. 1 2008b using a dust opacity 
of 0.01 cm 2 g~ 4 . lOssenkopf & Henning Il994h from the idealised 
GC (§ O- The b lue dotted line i n Figu re [3] shows the expected 
1mm flux from the lD’Ercole et ali l2008 ) model, when extrapolat¬ 
ing the density and temperature profiles in § !2.2l out to a radius of 
250pc (the maximum radius of the model in their Figure 1). The ex¬ 
pected large reservoirs of gas should be readily detectable in young 
massive clusters in galaxies out to several tens of Mpc. 

Finally, we consider what affect variations in the gas distri¬ 
bution would make to the above calculations. The uniform volume 
filling factor assumed in < 12. II represents the lowest possible peak 
extinction and the lowest possible gas and dust surface brightness. 
Due to mass conservation, the total amount of gas within the vol¬ 
ume considered must remain the same, however it is distributed. 
Therefore, a higher degree of gas clumping (equivalent to a smaller 
gas volume filling factor) will tend to make any extinction features 
more pronounced. In other words, although the area subtended by 
the extinction will decrease, the actual value of the extinction wher¬ 
ever the gas is confined should increase. 

If the gas distribution is not uniform, the spatial distribution 
of dense clumps can affect the optical/IR extinction and luminosity 


5 We note that iDalessandro et ai] i20li) observe multiple stellar pop¬ 
ulations towards a cluster which they estimate had an initial mass of 
~ 2 x 1O 5 M0. If this turns out to be common, evidence of centrally- 
concentrated reddening/extinction may be expected towards clusters much 
less massive than 10 6 M©. 
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Figure 2. Predicted extinction in the V-band in magnitudes based on the gas profile in the lD’Ercole et al.1120081) model (see text for details). The left and right 
panel show the extinction for the inner 30 pc and inner 8 pc, respectively. The colour scale, denoted by the colour bar, accentuates the intermediate extinction 
regions at each spatial scale. In the left panel, the outer, intermediate and inner black contours show the radius at which the extinction equals 0.1, 1 and 10 
mag. In the right panel, the contours show extinctions of 0.5, 1, 5, 10, 50 and 100 mag. 


profile. For example, although highly improbable, if all the dense 
clumps happened to lie on the far (near) side of the cluster along our 
line of sight, then they would sit behind (in front of) the majority 
of stars. In which case they would block less (more) of the cluster 
light than if they sat on the opposite side of the cluster along our line 
of sight. This would result in a less (more) pronounced dip in the 
optical/IR luminosity profile and less (more) pronounced optical/IR 
extinction features. 

Therefore, combining these effects, if the optical/IR obser¬ 
vations can spatially resolve the size of the clumps, and they 
do not all lie behind the stars, the clumps will be detectable as 
very pronounced (much higher extinction than calculated in Fig- 
ureQ}, patchy extinction features. If the observations cannot resolve 
the typical size scale of the clumps, the expected extinction and 
gas/dust emission becomes more complicated to determine. The 
optical/IR luminosity profile will then depend on both the size dis¬ 
tribution and location of the clumps. Although certainly interesting, 
these more detailed calculations are beyond the scope of the current 
paper. 

Unlike the case for optical/IR extinction, the emission from 
the gas and dust at mm wavelengths will be the same whether the 
clumps lie at the near or far side of the cluster. If the clumps are 
resolved, the higher gas volume density will lead to a much higher 
gas and dust surface brightness (i.e. brightness temperatures), mak¬ 
ing the gas and dust much easier to detect than shown in Figure [3] 
If the clumps are not resolved, the observed brightness temperature 
of the dust and gas emission will be reduced by approximately the 
“bearr@ filling factor” - the total area subtended by all the clumps 
within the beam, divided by the beam area. 

In summary, if observations can resolve the typical size scale 
of any gas substructures, the assumption of uniform volume fill¬ 
ing factor represents a conservative estimate of the expected peak 


The ‘beam’ is effectively the resolution of the observations. 



Distance [kpc] 

Figure 3. Expected flux density at a wavelength of 1mm for gas of a given 
mass, (dust) temperature and metallicity as a functio n of distance. The blue 
dotted line shows the predicted flux based on the D’Ercole et alj fe008l) 
model (see text for details). The horizontal dashed line shows the approxi¬ 
mate ALMA lOcr sensitivity limit at this wavelength for a one hour observa¬ 
tion. The vertical dotted lines show the distances to some nearby galaxies. 
Large gas reservoirs should be readily detectable towards young massive 
clusters in many nearby systems. 

extinction and minimum expected gas and dust brightness temper¬ 
ature. 


4 DISCUSSION 

Given current instrumental limitations, it would be extremely dif¬ 
ficult to detect the presence of gas reservoirs from the luminos¬ 
ity profile or mm continuum emission in young globular clusters 
at high redshift. However, these features should be readily de¬ 
tectable in clusters of similar stellar mass and density with ages 
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<100Myr in nearby galaxies with metallicities >0.1 Zq. Neither 
dips in the optical/infrared luminosity profile nor mm continuum 
emission have been detected towards any such young massive 
clusters (YMCs) fe.g. iBastian & Strader! 120 14l : ICabrera-Ziri et al.l 
l2015t) . We now consider potential explanations for this. 

A fundamental assumption of the preceding analysis is that 
the dust properties (grain size distribution, opacities etc) are simi¬ 
lar to those observed in the Milky Way and nearby galaxies (e.g. the 
large and small Magellanic clouds). If the dust properties were dif¬ 
ferent, this would be detectable through a combination of radio and 
sub-mm wavelength observations with optical/infrared extinction 
observations. If the grains were predominantly much smaller, one 
would expect to see copious PAH emission at 8/j.m in an environ¬ 
ment with such high Lyman-Werner luminosity. To our knowledge, 
bright PAH emission has not been observed towards the centre of 
YMCs. Extinction observations are typically not sensitive to dust 
grains much larger than the observing wavelength, so it may be pos¬ 
sible to hide the dust if the grains are all much larger than ~mm. 
However, in the absence of any plausible physical mechanism to 
cause all the dust grains to be >mm-sized, this seems highly im¬ 
probable. We dismiss anomalous grain size distribution variations 
as an explanation for the lack of dust detected towards young mas¬ 
sive clusters. 

A second assumption is that the gas-to-dust ratio scales in a 
similar way with metallicity as observed in th e Milky Way, Mag - 
ellanic clouds and other star-forming galaxies iFisher et all20141 ). 
If the gas-to-dust ratio were anomalously large, it may be possi¬ 
ble for clusters to contain large gas reservoirs without any signs of 
extinction or emission. However, such large gas reservoirs may be 
detectable in other ways. Once the gas becomes self-gravitating, 
cools to become molecular and begins forming stars it should be 
extremely bright in molecular gas tracer s such as C O and easily 
detectable with mm interferometers (see lBolatto et al.ll2013l) . The 
presence of gas with mass of order the stellar mass may also af¬ 
fect the dynamics of the stars, leading to larger-than-expected dy¬ 
namical mass to light ratios. No molecular line emission or anoma- 
lou s mass to light ratios have been observed towards YMCs (e.g. 
ICabrera-Ziri et alj|2015i) . Again, in the absence of any plausible 
physical mechanism to cause this, we dismiss anomalous gas-to- 
dust variations as an explanation for the lack of gas and dust de¬ 
tected towards young massive clusters. 

An alternative explanation is that the gas accumulation and 
star formation may happen in very short time intervals. This would 
mean that the chance of catching any individual YMC between 30 
and 100 Myr with detectable amounts of gas would be small. In or¬ 
der for this to happen, a mechanism is required to stop the mass 
being expelled from AGB stars for a long period of time. The re¬ 
quired mass reservoir from AGB stars with the correct abundances 
would then need to be released and migrate to the cluster centre on 
a very short timescale, as soon as the required reservoir of pristine 
material is being accreted. Such a scenario seems highly implau¬ 
sible. Nevertheless, if we take the extreme ‘burst’ scenario, where 
the gas is visible for 1 Myr over a 70 Myr period, we may only ex¬ 
pect to see 1/70 YMCs in the appropriate age range with signs of 
a large gas reservoir. This may provide an explanation for why ob¬ 
servations of individual YMCs of sufficient mass do not_show any 
evidence for reddening or extinction (e.g. McCradv et alj[2005l) . 
However, a direct prediction of this scenario is that observations 
of large samples of YMCs should recover some sources with ex¬ 
treme extinction/reddening of the luminosity profile. In addition, 
a consequence of such a large second star formation burst, is that 
some clusters in this mass and age range should show evidence for 


signs of extende d star formation histories. Recent studies find no 
such evidence IBastian et al.ll2013l : ICabrera-Ziri et alj2014l) . 


5 CONCLUSIONS 

If clusters of mass >10 6 Mg and radii of a few pc commonly form 
a second generation of stars at ages between 30 and 100 Myr from 
a large gas reservoir that builds up in the cluster centre, evidence 
for this gas would be straight forward to observe. The fact that no 
such observational evidence exists calls into question whether such 
a mechanism happens regularly for YMCs in galaxies within a few 
tens of Mpc. Future spatially-resolved optical/infrared and mm- 
continuum observations of large samples of YMCs are required to 
determine the fraction of such clusters with evidence of large gas 
reservoirs and/or recent/ongoing star formation. Such observations 
would be able to quantitatively test the predictions of GC formation 
scenarios that require GCs to accrete large gas reservoirs in order 
to form the second generation of stars. 
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